Micro-mechanical oscillator ground state cooling via intracavity optical atomic excitations 
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We predict ground state cooling of a micro-mechanical oscillator, i.e. a vibrating end-mirror of an optical 
cavity, by resonant coupling of mirror vibrations to a narrow internal optical transition of an ensemble of two 
level systems. The particles represented by a collective mesoscopic spin model implement, together with the 
cavity, an efficient, frequency tailorable zero temperature loss channel which can be turned to a gain channel of 
pump. As opposed to the case of resolved-sideband cavity cooling requiring a small cavity linewidth, one can 
work here with low finesses and very small cavity volumes to enhance the light mirror and light atom coupling. 
The tailored loss and gain channels provide for efficient removal of vibrational quanta and suppress reheating. 
In a simple physical picture of sideband cooling, the atoms shape the cavity profile to enhance/inhibit scattering 
into higher/lower energy sidebands. The method should be applicable to other cavity based cooling schemes for 
atomic and molecular gases as for molecular ensembles coupled to stripline cavities. 
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Cavity cooling of micro-mechanical oscillators has seen 
tremendous experimental progress IHHjlUHIllSSlMl in the 
past few years, also stimulated by theoretical analysis showing 
that the quantum ground state of such oscillators is achievable 

[l^ [Till . In essence, through cavity mediated interaction 
the mechanical vibrations at kHz-MHz frequencies are con- 
verted to strongly damped sidebands of light modes at opti- 
cal frequencies, that contain essentially no thermal excitations 
(thermal photons). Standard sideband cooling is effective only 
when the linewidth of the optical resonator k is smaller than 
the mirror vibration frequency ct;,„ and one has a sufficient op- 
tomechanical coupling. For a given mirror reflectance a small 
enough k can be reached using longer cavities, but strong op- 
tomechanical coupling then requires a very large intracavity 
field and input power. 

Here we propose the alternative use of atomic, molecular or 
solid state electronic transitions as extremely localized narrow 
bandwidth oscillators, for the cooling process. These are cou- 
pled to the vibrational modes via a cavity light mode. Strong 
mirror-light and light-atom coupling can be achieved while 
still maintaining a narrow and tailorable resonance behavior 
As the simplest generic example we will consider ensembles 
of two-level atoms with a narrow transition suitably detuned 
from the cavity resonance. At optical frequencies thermal ex- 
citations of the atoms can be neglected and the atoms consti- 
tute an effective zero temperature reservoir, where energy and 
entropy can be dissipated efficiently. Note that using optical 
resonances to cool solids has already been shown successful 
in optical fibers 1,1^ ^. 

Cavity assisted sideband cooling is based on unequal 
scattering of pump photons into higher/lower energy (anti- 
Stokes/Stokes) photons; the cavity modifies the density of 
optical modes around the laser frequency such that a higher 
density at the anti-Stokes frequency leads to a higher scat- 
tering rate into the more energetic sideband. Denoting the 
corresponding rates with Aas and Ag, an effective extrac- 
tion rate of vibrational quanta T = Aas — Ag can be added 
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FIG. I: (a) Typical resolved-sideband cooling of a micro-mirror in 
the good-cavity regime is achieved via cavity enhancement of the 
blue (anti-Stokes) sideband. Laser driving is at frequency coi while 
sidebands are scattered at cos and ujas- The dashed black line shows 
the cavity profile around the ujc cavity resonance, (b) Cooling via 
two-level atoms in the bad-cavity regime. The empty cavity response 
around coc (dashed black line) is modified by the presence of the 
atoms of frequency splitting uji which leads to a dip at ujs (green 
line), thus effectively inhibiting Stokes scattering (heating), (c) Cool- 
ing via inverted two-level atoms with splitting uj2. The peak at uJas 
(green line) leads to enhancement of the anti-Stokes scattering (cool- 
ing). 



to the intrinsic decay rate of the mechanical oscillator 7^ 
to yield j„i = F + 7,„ ^ 7^. For an oscillator of fre- 
quency (jj„i, this produces cooling to an effective tempera- 
ture T ~ /ia;„i/ [fcs ln(l + 1/n)] (ks =Boltzmann's con- 
stant) where fi = •jmn/'jm + rires- Therefore cooling leads 
to a reduction of the initial average thermal vibration quantum 
number 77 by a factor 7„i/7m, but one has also an additional 
residual occupancy rf''^'' ~ Ag/^m that is due to the residual 
Stokes scattering rate. Hence one best uses a high-finesse op- 
tical cavity (with k < ujm, as illustrated in Fig.[T^) resonant 
with the blue anti-Stokes sideband [loll leading to cooling 
ton'^'^'* ~ {n/2ujmf. 
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The novel mechanism that we propose in this Letter can be 
well understood from the picture of cavity cooling shown in 
Fig.[TH and is illustrated in Fig.[T}5,c. The added ensembles 
of two-level atoms, one prepared in the ground state (type 1), 
and the other one inverted (type 2), modify the cavity suscep- 
tibility, and consequently change the cavity induced scatter- 
ing to the sidebands. Even in the limit of a cavity bandwidth 
much larger than the mirror vibration frequency the atoms 
still change the response function to inhibit Stokes scattering 
and enhance scattering into the anti-Stokes sideband. As we 
will see later, the first mechanism is much more efficient than 
the second in that, while it is characterized by slower cooling 
rates, it ensures a low n^'^'^ <^ 1, that allows for ground state 
cooling. The difference comes from the fact that while the 
cavity filled with ground state atoms forms an effective zero- 
temperature bath, the noise inherent in any amplification pro- 
cess implies an effective non-zero temperature near the gain 
maximum of the inverted atoms. Cooling in the second case 
is then bounded by the effective bath temperature. 

Hamiltonian - A Fabry-Perot optical cavity is driven by a 
laser of frequency loi near one of the cavity resonances lUc 
with detuning Ao ~ lOc ~ uJi- One of the cavity mirrors is 
allowed to vibrate at the mechanical frequency ti>,„ and with 
an effective vibrating mass m. Two clouds of iVi 2 two-level 
atoms (TLA) with level splittings ^1,2 (detuned from the laser 
by Ai_2 = ti'1,2 ~ '-^i) are positioned inside the cavity mode. 
The Hamiltonian of the system is {h — 1) 

H = Ho + Hi + Hms, (la) 
Ho = Uca^ a + [p^ + q^) + uiSiz + UJ2S2Z (lb) 

Hi = gi {Si+a + Si-a^) + g2 {S2+a + S2-a^) (Ic) 
-Goa^aq + i£ (a^e"*'^'* - ae*""*) , 

where Hq is the free part, Hj is the interaction part while 
Hdis is the dissipative part which will not be explicitly shown 
here but included in the Langevin equations. The free part 
includes the cavity mode, mirror and atom free energy where 
a is the annihilation field operator with [a, at] = 1, q and 
p are the mirror quadratures with [q,p] — i and the atoms 
are described by a set of collective spin operators satisfy- 
ing [5i,2+,S'i,2-] = 5*1^2^ and [Siaz,Sia±] = ±S'i,2±. 
The interaction contains the atom-field Tavis-Cummings term 
and the radiation pressure field-mirror coupling where 51 2 
is the single atom-field coupling strength while Go = 
(lOc/L) y^h/mujm (L is the cavity length) is the single photon- 
mirror optomechanical coupling. The laser driving shows up 
in Hj as a displacement term of amplitude £ (directly con- 
nected to the input laser power). Not included in the Hamil- 
tonian is a driving term that leads to an effective population 
inversion in type 2 atoms. 

Linearization - For simplicity we perform a transformation 
from the spin algebra to a harmonic oscillator algebra for the 
atomic system via ci = Si-/^/Ni and C2 = ^2+ /VN2, such 
that: ujiSiz ~ uji{-Ni/2 + c{ci) and LJ2S2Z ^ ^2(^^2/2 - 
c|c2) and we have [ci,cj] = (5y. The effect of H^is in 
Eq. ([T]i can now be added in a Langevin equation approach 
via the quantum noises ai„, pin, ci.i„ and C2.m with the fol- 



lowing nontrivial correlations: (am(i)ai„(i')) = 2K5{t — t'), 

{p^n{t),P^nit')) = 7m(2n + - {ci MM ,^nit' )) = 

2ji5{t ~ t') and (c2,,„(t)4,,„(t')) - 272(5(t - t') Q. The 
cavity and atom reservoirs are at zero-temperature and are 
fully characterized by the decay rates k, 71 and 72 while the 
thermal mirror reservoir is described by the damping 7„i and 
thermal occupancy n. Notice that 71 is the actual decay rate 
of atoms of type 1 to their real ground state; 72 is instead the 
effective optical pumping rate to the excited state and can be 
considered as an effective "decay" rate. The other levels in- 
volved in the optical pumping process are irrelevant for cool- 
ing and are neglected in the present two-level description. 

As in 1 14!], one can linearize the equations of motion for the 
atom-field-mirror system around steady state values to derive 
a set of quantum linearized Langevin equations (QLLEs) 



q = uj^p, (2a) 

p = -uj„iq - 7,„_p + G (a + a^) +pin, (2b) 

d = — {k + iAf) a + iCq — iGiCi — iG2c\ + a,i„, (2c) 

ci = - (71 + zAi) ci - zGia + ci,i„, (2d) 

C2 = - (72 - i^2) C2 - iG2a^ + C2,m. (2e) 



where all the couplings are collectively enhanced, Gi,2 = 
.91,2 ^1,2 and G = Goa, with a the steady state intracav- 
ity field amplitude (which can always be assumed real). In 
Eqs. ^ a, q, p and ci^2 operators denote steady state fluc- 
tuations, while Ay = Ao — G"^ jijJm is the effective cavity 
detuning. 

In the following it will be useful to define 

l/[K + i(A/-tj)], ei(tj) = 1/ [71 + i (Ai - tj)] and 
£2 (w) ~ 1/ [72 — * (A2 + w)] as the mirror, field and atom 
bare susceptibilities. Fourier transforming Eqs. (O one can 



derive the following set of coupled equations 

[em(t^)]"^g = G(a + aT) +]3j„, (3a) 

[e/ {yi)Y^ a = iGq + d^, (3b) 

where the effective field input noise is din = o-in — 



i[Giei{ujt)ci,in + G2f-2[—^)c2 in]- Th^ fisld response func- 
tion is modified by the atoms according to 

[If {uj)]-' = [ef {uj)]-' + G?ei {uj) - GY2 (-^) • (4) 

Physical picture. Cavity response - In the absence of atoms 
(setting G1.2 = 0) Eqs. ([3]) describe a typical optomechani- 
cal setup, which in the good-cavity limit k <C ojm (resolved- 
sideband limit) and A/ = leads to optimal cooling of 
the mirror state. This effect can be simply understood in 
terms of preferential scattering of sidebands at frequencies 
ijJas,s ~ ^i^^m off the mirror, as illustrated in Fig.[T^, where 
the cavity response (i.e., |e/ (cj)| is plotted around Wc) to the 
mechanical sidebands is depicted. The choice of lOc ~ oJas 
allows the enhancement of the cooling sideband. When atoms 
are placed inside the cavity, their effect on the cavity response 
function leads to a modification of the cavity mode structure 



3 



around the resonance ujc quantified by Eq. (|4|i. This is illus- 
trated in Fig.[TJ) for the bad-cavity limit where k ^ w^, and 
resonant cavity driving loi = lUc- The type 1 atoms, placed 
at Ai = — cj„i induce a dip in the cavity profile at lOs while 
the gain medium increases the cavity response at Was- What 
results is an inhibition of scattering processes into the Stokes 
sideband and enhancement of scattering into the cooling side- 
band. 

Based on the mode structure tailoring picture, we can now 
make some quantitative remarks on the efficiency of cooling 
that is expected from such a scheme. With the usual defini- 
tion of the atom cooperativity Ci^2 = G\ 2/ (^71, 2), one can 
see that the dip value at ujs scales down as (1 + Ci) ^ while 
the peak at LOas has a maximum that scales up (1 — C2) 
The widths of the dip/peak are 71 = 71 (1 + Ci) and 72 = 
72 (1 — C2), respectively representing the enhanced/reduced 
light -induced atomic line widths. In the limit where the scat- 
tered sidebands fit inside the dip/peak, i.e. the resulting cool- 
ing rate is smaller than 71, 72, one then expects an inhibition 
of the Stokes peak (enhancement of the anti-Stokes peak) by a 
factor of the order of(l + Ci)~^ ((1 — (^2)^^, respectively). 

Cooling - For a proper derivation of cooling rates and resid- 
ual occupancy we start with Eqs. Q. Notice the similarity 
with typical cavity-assisted cooling of a mirror, which can 
be obtained as a limiting case when Gi_2 = for which 
o-in = a-in and If {uj) = £/ {ij->) ifiol [Till . The mechanical 
oscillator of Eq. ([3^) is driven by a thermal zero-average 
Langevin force pi„, and by a field-atom induced Langevin 



force F = G [ain + a{„ j . Following the approach of iK 
one can find the cooling rate induced by the field-atom system 
and the mirror occupancy directly from the spectrum of the 
force F, (F (w) F (w')) = Sp (w) 5 {uj + uj'). From Eq. ©) 
one derives 
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+ |e-/(-c.)|S2G^|e2Hr}. (5) 



The cooling and heating rates (anti-Stokes and Stokes, respec- 
tively) are computed as Aas.s = Sp (ii^m) /2. We shall con- 
sider from now A/ = 0, i.e., the cavity resonant with the laser, 
so that there is no cooling in the absence of atoms. 

Cooling with ground state atoms - We restrict our discus- 
sion in the following to a cavity filled with atoms of type 1 
(setting G2 = 0) and derive the effective cooling rate and 
residual occupancy under resonance condition Ai = ~ujm- 
First we stress that the effective temperature of the cavity is 
zero: {a\^{ijj)ain{ijj')) = 0. We identify the optimal regime 
as the bad-cavity regime k ^ w™, where the effect of the 
atoms on the cavity profile is to efficiently inhibit Stokes 
scattering for a bandwidth 71 around w,,. We ask that the 
bandwidth 71 be larger than the scattered sideband halfwidth, 
which is given by the total mechanical damping rate Fi + 7,„ . 
In this limit and with the requirement 71 ^ cj„i,. one can 
approximate \ei{urn)\ ^ 1/ (2u;„i), |ei(-Wm)| = I/71. 
\ef{u:ra)\ ^ 1/k and |e/(-a;„)| ~ + The 

scattering rates can be derived from Eq. (|5]) as Aas — G"^ / n 
and As ~ G^/ [k (1 + Gi)] and consequently the cooling rate 
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FIG. 2: Relative cooling rate 7^/ 7„i and residual occupancy n^es 
(in the insets) as a function of the normalized parameters Gi /um 
and G2/ijJm for G = uim- The red line corresponds to analytical 
expressions while the black line shows the exact results (see text), 
a) Cooling with ground state atoms and G2 = 0. b) Cooling with 
inverted atoms and Gi = 0. The other parameters are L0m/2n — 10 

MHz, 7,„ = 10"^ X UJm, K, = 100 X UJ„„ 71 = LJm/100, 72 = CJm 

and n = 100. 



G^ Gi 
~1 + Gi 



(6) 
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The residual occupancy can be computed as n'l 
As/j,n — As/Ti Gi^ (in the limit of large coopera- 
tivity). This result replaces the resolved-sideband limit final 
occupancy rf'^^ = (K/2a;,„)^ of the purely optomechanical 
case with no atoms. The advantage of adding ground state 
atoms within the cavity is obvious from the scaling of rfi"^^ 
with the controllable and in principle unbounded parameter 
Ni. In fact, in the strong resolved-sideband limit k ^ lo^, 
the cavity is far-off resonance, and one has to use large input 
power for optimal cooling, while with atoms one can make 
smaller and smaller by increasing A^i . Illustrations of the 
dependence of ^m/lm and on the normalized coupling 
strengths Gi/wm when G2 = is shown in Fig.|2^ where red 
lines correspond to the exact numerical solution of Eqs. (|2]) 
in steady state to derive the 4-mode covariance matrix, while 
black lines show the corresponding previously derived analyt- 
ical expressions. The set of parameters we have made use of 
is cj„i/27r = 10 MHz, 7,„ = 10^^ x ujm, n ~ 100 x uj^, 
71 = uj„JlOO and n = 100. 

The condition Fi + 7„i < 72 puts a lower bound on 
Gi which is given by G. Moreover, with the increase of 



Gi, the field response at ujas (|e/ ) diminishes as 1 + 

(Gf /2Ktj„i)^ from the value of 1/k^ assumed above, which 
in consequence amounts to the limitation that Gi < y/KOJm. 
The cooling rate is independent of Gi as long as Gi > G, 
however the residual occupancy is limited by 71^^ ~ Gi^ > 
7i / uJm, which is the justification for working with small 71 . 

Cooling with inverted atoms - We shift our discussion now 
to the opposite case of a cavity filled with atoms of type 2, 
at A2 = uj„i and set Gi = 0. The correlations of the input 
noises are now (aj„(w)am(w')) = 2j2G2\e2{uj)\'^S{iL! + lo') 
and (ai„(a;)aj„(a;')) = 2k5{ijJ + oj') which show that the 
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cavity is at a nonzero effective temperature. This results in 
a frequency-dependent thermal cavity occupancy n (uj) = 
72G||e2(w)pK"^(l - 72G'||e2(a;)P/K)"^ which evaluated 
at ujjn leads to n{uj„i) ~ C2I (1 — C2). This will set a 
lower bound for the achievable mechanical occupancy be- 
cause the oscillator cannot be cooled to less than the occu- 
pancy of the effective cooling reservoir. In the bad cavity limit 
the atoms produce a peak of bandwidth 72 around to as- We 
ask again for the condition that the anti-Stokes sideband fits 
inside the peak, i.e. r2 + 7m < 72- In this limit and with 
the requirement 72 <C i/J^-^, one can approximate |e2 — 
l/(2t^m), |e2(-cc;,„)| = 1/72, |e/(t^m)| ^ 1/ [(1 - C2) ^t] 
tJm)! ~ The scattering rates are Aas — 



and 



and 



and the cooling rate is 



GVk 1 + G2/(1-G2)^ 



G^ G2 

K 1-G2 



(7) 



The residual occupancy in the limit of G2 ^ 1 diverges as 
^res ^ G2/ (1 — G2), which is the occupancy of the effec- 
tive reservoir at the Stokes sideband. However, G2 cannot 
be increased arbitrarily close to 1 since there is a constraint 
r2+7m <72 which leads to G^/ (72 k) < (1 - G2)^. Graph- 
ical illustration of the increase of the normalized cooling rate 
1ml Im as a function of G2 jiOm at fixed G\ = is provided in 
Fig-Eb- The analytical result is valid up to values of G2 where 
7m is of the order of 72 . The inset shows the agreement be- 
tween the analytical and numerical results for ril^^ when G2 
is large such that G2 ^ 1. 

Discussion - The use of ground state atoms for cooling 
presents the advantage of a low residual occupancy propor- 
tional to Gf ^. The downside is that cooling is slow at a 
rate that saturates at G^/k in the limit Gi ^ 1. The use 
of inverted atoms considerably speeds up the cooling pro- 
cess that takes place at a rate r2 ^ G'^/k for G2 close 
to unity, at the price, however, of a high rires- One would 
ideally combine the two methods for optimal cooling. How- 
ever, such a mixed scenario is not particularly advantageous. 
In fact, in the limit of 71,2 <C cj„i, one obtains a combined 



cooling rate that is practically the one induced by the anti- 
Stokes enhancement ri2 ~ r2, but the residual occupancy 
^r-es ^ n^:^<'(\ — C2) + n^-i^ ~ 712'^* is still large, reflecting 
the negative effect of the high effective temperature induced 
by the inverted atoms. 

A few remarks on the system stability have to be made at 
this point. The analysis can be made by applying the Routh- 
Hurwitz criterion II15I1 on Eqs. A cavity with a moving 
mirror driven at resonance is always stable and its stability is 
not perturbed by the addition of ground state atoms at Ai = 
—ijJra- The inverted atoms, however, do impose a limitation 
on the system's stability, which can be simply cast as G2 < 1 
and trivially fulfilled by the proper choice of 

The question of the validity of bosonisation also has to be 
discussed. For the ensemble of two level atoms to resemble 
an harmonic oscillator one has to ask that the single atom 
excitation induced by the cavity field is small. This means 
.91.2/ (^m + 71.2) < < 1' which, as detailed in iH] 
can be fulfilled for a dipole forbidden transition or for atomic 
vapor cell much larger than the cavity mode. 

Conclusions - We have shown a new mechanism for mir- 
ror ground state cooling where atoms are used to dissipate the 
thermal energy of the mechanical oscillator. In a simple pic- 
ture this can be seen as cavity-assisted cooling with a cavity 
profile tailored via use of atoms at the Stokes and anti-Stokes 
sidebands. The effect is to inhibit scattering processes that 
lead to heating while enhancing the cooling process. We no- 
tice that there are alternative ways of cooling through cavity 
spectrum tailoring. A recent example is given in [16] where 
tailoring is obtained when the micro-mechanical resonator di- 
rectly modulates the cavity bandwidth. More generally speak- 
ing, appropriate spectrum tailoring can be always be obtained 
through quantum interference |17i], as for example proposed 
and demonstrated for cooling trapped ions ifisll . 

Acknowledgements - We are grateful to G. Milburn and 
M. Wallquist for helpful discussions. We acknowledge finan- 
cial support from the European Commission (FP6 Integrated 
Project QAP, and FET-Open project MINOS) and Euroquam 
Austrian Science Fund project II 19 N16 CMMC. 



[1] C. H. Metzger and K. Karrai, Nature (London), 432, 1002 
(2004). 

[2] O. Arcizet, P.-F. Cohadon, T. Briant, M. Pinard, and A. Heid- 

mann, Nature (London) 444, 71 (2006). 
[3] D. Kleckner and D. Bouwmeester, Nature (London) 444, 75 

(2006). 

[4] T. Corbitt, Y. Chen, E. Innerhofer, H. MuUer-Ebhardt, D. Ot- 

taway, H. Rehbein, D. Sigg, S. Whitcomb, C. Wipf, and N. 

Mavalvala, Phys. Rev. Lett. 98, 150802 (2007). 
[5] J. D. Thompson, B. M. Zwickl, A. M. Jayich, F. Marquardt, S. 

M. Girvin, and J. G. E. Harris, Nature (London) 452, 72 (2008). 
[6] J. D. Teufel, J. W. Harlow, C. A. Regal, and K. W. Lehnert, 

Phys. Rev. Lett. 101, 197203 (2008). 
[7] A. Schliesser, O. Arcizet, R. Riviere, T. J. Kippenberg, 

larXiv:0901.1456v l 



[8] S. Groeblacher, J.B. Hertzberg, M.R. Vanner, S. Gigan, K.C. 

Schwab, M. Aspelmeyer, arXiv:0901.1801 A 
[9] I. Wilson-Rae, N. Nooshi, W. Zwerger, and T. J. Kippenberg, 

Phys. Rev. Lett. 99, 093901 (2007) 
[10] F. Marquardt, J.P Chen, A.A. Clerk, and S. M. Girvin Phys. 

Rev. Lett. 99, 093902 (2007). 
[11] C. Genes, D. Vitali, P. Tombesi, S. Gigan, and M. Aspelmeyer, 

Phys. Rev. A 77, 033804 (2008) 
[12] A. Rayner, M.E.J. Friese, A.G. Truscott, N. Heckenberg, and 

H. Rubinsztein-Dunlop, J. Mod. Opt. 48, 103(2001). 
[13] C. W. Gardiner and P. ZoUer, Quantum Noise, (Springer, Berlin, 

2000). 

[14] C. Genes, D. Vitali, and P Tombesi, Phys. Rev. A 77, 050307 
(2008). 

[15] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and 



5 



Products, Academic Press, Orlando, pag. 1119, (1980). 
[16] F. Elste, S.M. Girvin, and A.A. Clerk, Phys. Rev. Lett. 102, 

207209 (2009). 
[17] D.A. Rodrigues, Phys. Rev. Lett. 102, 067202 (2009). 
[18] G. Morigi, Phys. Rev. A 67, 033402 (2003); G. Morigi, J. Es- 



chner, and C.H. Keitel, Phys. Rev. Lett. 85, 4458 (2000); C.F. 
Roos, D. Leibfried, A. Mundt, F. Schmidt-Kaler, J. Eschner, 
and R. Blatt, Phys. Rev. Lett. 85, 5547 (2000). 



